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Abstract 

We review aspects of classical and quantum mechanics of many anyons confined in an oscil- 
lator potential. The quantum mechanics of many anyons is complicated due to the occurrence 
of multivalued wavefunctions. Nevertheless there exists, for arbitrary number of anyons, a 
subset of exact solutions which may be interpreted as the breathing modes or equivalently 
collective modes of the full system. Choosing the three-anyon system as an example, we also 
discuss the anatomy of the so called "missing" states which are in fact known numerically and 
are set apart from the known exact states by their nonlinear dependence on the statistical 
parameter in the spectrum. 

Though classically the equations of motion remains unchanged in the presence of the statis- 
tical interaction, the system is non-integrable because the configuration space is now multiply 
connected. In fact we show that even though the number of constants of motion is the same 
as the number of degrees of freedom the system is in general not integrable via action-angle 
variables. This is probably the first known example of a many body pseudo-integrable system. 
We discuss the classification of the orbits and the symmetry reduction due to the interaction. 
We also sketch the application of periodic orbit theory (POT) to many anyon systems and 
show the presence of eigenvalues that are potentially non-linear as a function of the statistical 
parameter. Finally we perform the semiclassical analysis of the ground state by minimizing 
the Hamiltonian with fixed angular momentum and further minimization over the quantized 
values of the angular momentum. 
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1 Introduction 



In 1977, Leinaas and Myrheim^ showed that in two space dimensions it is possible to have 
particles obeying intermediate statistics different from the well known Bose and Fermi statistics. 
Wilczekj2] later coined the name anyons for particles obeying these peculiar statistics. In the 
last two decades, a lot of work has been done in this field. Many of these developments have 
been brought up to date and nicely summarized in the books by Lerda[Sl and KhareQ as also 
in the collection of articles on anyon superconductivity 0. One aspect that has not been covered 
in detail in these books and other reviews^ concerns the question of integrability properties of 
a system of many anyons and its classification. This review is intended to focus on these aspects 
which have not been discussed in detail earlier or have been largely left out of earlier attempts. 

Here we focus on the non-relativistic classical and quantum mechanics of particles in two 
spatial dimensions which obey fractional statistics ^H]. Generically such particles are referred 
to as anyons. Anyons as a research area is now over twenty years old. These systems have emerged 
as being interesting in their own right from the point of view of mathematical physics at both 
the classical and the quantum level. These systems also constitute an example of inequivalent 
quantizations due to non trivial fundamental group of the configuration space 7\. It is well 
known that the existence of fractional statistics is intimately connected to having multivalued 
wavefunctions which naturally occur in quantum mechanics on multiply connected spaces jHj. 
The classification of multivalued wave functions is provided by one dimensional representations 
of the fundamental group of the multiply connected configuration space. The configuration space 

is the d-dimensional Euclidean space given by R2 with all the diagonal points, A, removed. 
The fundamental group of turns out to be 

TTiiQ^) = Sn if d > 3 

and 

7ri(Q^) = 5,v if d = 2 
where Bj\f is the braid group of N objects which contains the permutation group Sn as a subgroup. 
This immediately brings out the difference between two and three dimensions. Although there 
are several known examples of a kinematic classification [7j (in terms of the representations of the 
fundamental group), very little is known about the dynamics on such spaces. The motivation for 
studying such systems is two fold: they are a good and physically relevant example of quantum 
mechanics on multiply connected spaces. It is physically relevant because anyons were proposed 
as candidates for explaining fractional quantum Hall effect [HI and their possible connection to 
high Tc superconductivity [TOj. It is a good example because at least some exact solutions to 
the energy eigenvalue problem for N anyons in some external confining potential are known — an 
exception in many-body quantum mechanics with non-separable Hamiltonians. 

The quantum spectrum of anyons has been analyzed by many authors. Broadly speaking, 
several authors have pointed out the existence of a subset of exact solutions TTI whereas in the 
numerical analysis of the low lying spectrum of three and four anyons, many non-trivial features 
of the spectrum have been shown to exist jl2j. The quantum mechanical spectrum (with harmonic 
potential added) shows two distinct qualitative features: a) eigenvalues which depend linearly on 
the statistical parameter, a, all of which are exactly known and b) eigenvalues which depend non- 
linearly on a and none of which is exactly known|12j. This is also borne out by a perturbative 
analysis of the three-anyon ground state I14j. Another interesting feature of the spectrum is 
the additional degeneracy around the semion point which is related to the reflection svmmetrvjlSj. 
These features may be traced to the properties of partial separability and pseudo-integrability 
manifested at the classical level jl6j . 
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Partial separability of the Hamiltonian may be exhibited explicitly in terms of two collective 
degrees of freedom and the remaining "relative" degrees of freedom. The total Hamiltonian, after 
removing the center of mass, may be written as a sum of Hi{coUective) + H2{collective, relative). 
The commutator of Hi and H2 turns out to be proportional to H2 JH]- This implies that the 
subspace of eigenstates of the full H on which H2 vanishes are exact eigenstates of Hi and these 
give all the exactly known eigenvalues and eigenstates. 

Pseudo-integrability, a concept first introduced by Richens and Berry jJT], may be described by 
exhibiting 2N constants of motion in involution, which fail to lead to integrability via usual action- 
angle variables. This was conjectured to be the reason for the level repulsion seen numerically in 
the non-linearly interpolating spectrum ^HI- The precise meaning of pseudo-integrability in this 
context, however, was not elaborated earlier. In particular, how this classical feature translates 
into there being only two good quantum numbers, was not analyzed. A further fallout of the 
pseudo-integrability property of the anyon system is its effect on the so called nearest neighbour 
spacing distribution in the energy spectrum. It is now known that the spectrum exhibits features 
of both regular and irregular spectra as a function of the statistical parameter |18j. 

Yet another method of approaching the quantum spectrum of many anyons is through the 
stationary phase approximation (SPA). In this, an approximation to the propagator G{E + ie) is 
developed using a suitable path integral representation ^21- In this approach the propagator is 
typically obtained as a sum over (families of) periodic trajectories in the classical phase space. This 
is entirely given in terms of classical quantities. One therefore expects to see directly the effects 
of non-trivial fundamental group of the phase space, pseudo-integrability and partial separability 
in a semiclassical framework. Since at the classical level, anyons with oscillator confinement, are 
locally identical to the oscillator system (but differ in the global topology of the phase space), the 
periodic orbits are known. Thus application of the periodic orbit theory (POT) to anyons, may 
be expected to be tractable. 

There is one clarification to be noted. At the quantum level, the statistical parameter, Og, 
enters via the stipulation of multivaluedness of the wave functions and is dimensionless. At 
the classical level its counterpart, ac, enters as the coefficient of a total derivative term in the 
Lagrangian and has the dimensions of h. These two must be related as ac = haq. In the SPA 
computation, Oc is held fixed with h going to zero. In the leading approximation, the spectrum 
will depend linearly on ac and thus also on Og. Alternatively, for comparison with the quantum 
spectrum to leading order, one will want to keep aq fixed implying that ac to go to zero with h. 
Viewed either way, one may not expect to see non-linear dependence on the statistical parameter 
at the level of the leading approximation. Nonetheless, if the non-linearly interpolating eigenvalues 
also have a linear piece, one should be able to see it in the leading SPA. 

This article is concerned with the discussion of the many points raised above. While many of 
these points have emerged from previously published papers, we attempt a coherent presentation 
of all these aspects with appropriate connections noted along with new material. We begin in 
section 2 with the basic formalism where we build the Hamiltonian of the anyons in two dimensions 
by requiring that their exchange statistics depend on a continuous parameter. In section 3, we 
review the quantum spectrum of many anyons. We first prove the existence of a subset of exact 
solutions. We then discuss, albeit in the case of three anyons, the so called missing states. 
Noting that eigenvalues themselves are independent of where the differential equation is solved, 
we expect the asymptotic form of the differential equation to determine the spectrum. The three 
anyon problem in the asymptotic regions reduces to solving a Confluent Hypergeometric Equation 
(CHE) in the region R < x < 00. This equation has a solution regular at x = and a solution 
which is irregular at x = 0. The regular solution leads to linearly interpolating energies while 
the irregular one may (R is strictly greater than zero) lead to nonlinear ly interpolating solutions. 
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Though we are unable to make definite predictions about the spectrum, we point out the inherent 
difficulty which involves the question "what is the quantization condition". 

In section 4, we analyze the many anyon system classically and show that the subset of exactly 
solvable solutions arise from the quantization of the collective coordinates (after the trivial centre 
of mass degrees are separated). This is made possible by the property of partial separability of 
the system for certain initial conditions and consistent with the equations of motion. In fact these 
solutions do not carry any information on the internal dynamics of anyons which may be frozen 
as far these solutions are concerned. These eigenvalues also goes as N'^ for N anyons. Thus the 
solutions are trivial dynamically. The real anyon dynamics therefore resides in the solutions which 
depend non-linearly on the statistical parameter. In the analysis of the quantum virial expansion 
for the equation of state of the anyon gas, it was shown that these trivial states conspire to 
cancel the divergent parts while the nonlinear solutions give the finite part which then defines 
the equation of state 20 . Thus it is of utmost importance to obtain these non-linear states to 
study the thermodynamic properties of the system. The low lying states have been obtained 
numerically ^21 and perturbatively |131 114j in the case of three and four anyons. In the case of 
large N anyonic systems, Thomas-Fermi approach has been used to obtain the properties of the 
ground and low lying excited states PT]. However to establish a virial expansion for the equation 
of state, by going beyond the third and the fourth virial coefficients, much more work on the 
nonlinear states is needed. 

In section 5, we discuss semiclassical approach via the periodic orbit theory, we consider 
the question of classical integrability- the equations of motion of course do not depend on the 
statistical parameter, however the configuration space of N- particles is R2 — A, where A denotes 
the set of all diagonal points which are removed. All the orbits which pass through these diagonal 
points are now modified. This leads us to conjecture that the system may be non-integrable a-la 
Richens and Berry jl7j . We demonstrate that the classical analysis is adequate for exhibiting 
quantum symmetries. In this process we sharpen and clarify the meaning of pseudo-integrability 
and role of the fundamental group of the phase space. The system is shown to be identical to the 
isotropic oscillator locally but not at the global topological level. We then classify the classical 
trajectories using symmetry transformations and explain how the non-trivial global topology 
reduces the symmetry group. Furthermore, we sketch the application of the periodic orbit theory 
to many anyons. In the process we show the presence of eigenvalues which are potentially non- 
linearly interpolating. As mentioned earlier such eigenvalues have been seen in the numerical 
spectrum for three and four anyons and have been conjectured to be present for general 
j22j . We obtain these systematically from POT. Further, both the classical modeling and the 
exact propagator for two anyons, indicate an ambiguity regarding possibility of closed orbits with 
half the basic period and its inclusion in the trace formula. We discuss this ambiguity in some 
detail. We show that it is possible to regularise the point-like statistical fiux and deduce the 
existence of half period trajectories. We look for their evidence in the exact propagator for two 
anyons and point out their ambiguous status. In Section 6, we present a detailed semiclassical 
analysis of the ground state of many anyons using a two step minimization. While this may seem 
unusual procedure to follow, we show that the method (even at the classical level) retains some 
of the features of the quantum spectrum. 

Finally the results are summarized in Section 7. The appendices contain many results to keep 
this review self contained. Appendix A contains a discussion of regularisation method for the 
so-called reflecting orbits which arise in the zero angular momentum sector. Appendix B contains 
details of OSp{4:N, R) symmetry. Appendix C contains details of the Symplectic diagonalisation 
procedure relevant for Section 6. 
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2 Basic Formalism 



In this section, we first develop the concept of anyons in two dimensions through multivalued 
wavefunctions. The exchange statistics of anyons then involves a continuous parameter, a — the 
so called statistical parameter^. It is the peculiarity of two dimensions that the quantization does 
not depend of the particular value of a, unlike in higher dimensions where we require a to be an 
integer. 

In the following, by "anyons" we mean a quantum mechanical system of N particles in two 
dimensions with wave functions which have a stipulated multi-valuedness to be specified below. 
To make this explicit let us denote a generic wave function as 'ip{fi, ...,fiy), where denotes the 
position vector of a particle. Let [Pijl-y denote the operation of taking the ith particle coordinate 
around the jth. coordinate along a closed path 7. The path 7 does not enclose any other particle 
coordinate and is taken in an anti-clockwise sense, say. Then let us stipulate that under such an 
operation ip acquires a phase namely, 



[Pij]^ip{fi,...rN) = exp{i2TTa)ip{fi, ...fiy)- 



(1) 



If a path 7 encloses other particle coordinates as well then such a path can be broken into a set 
of closed paths each of which encloses exactly one particle. Applying the stipulation above, one 
can compute the total phase, for such a path. If the sense of the path is reversed then a — > —a. 
Clearly the phase acquired depends only on the homotopy class of the path (i.e., it is the same 
for two paths 7 and 7' if 7 and 7' can be continuously deformed into each other). 



Let us introduce the complex notation for particle coordinates: zj 



+ iyj, Zj 



Zj, has the property that if Zj is taken around Zj, zfi; changes by 



Clearly Zij, where Zij = z 
exp(i27ra). This allows us to write any generic wave function satisfying eq.(^ as. 



'4'{Zi,Zi) 



n I 



a/2 



ljjiZi,Zi), 



(2) 



with the bracketed expression being a phase and now tl^^Zi, Zi) is a single valued function. 
Clearly, 



i<j 



which can be rewritten as, 



\ a/2 



Zij 



Vki^{zi,z.{) + Vkln n 



Zij 



a/2 



i'{Zi,Zi) 



(3) 



Vk'ip{zi,Zi) = 



i<j 



,\a/2 



Vki'izi^zi) + ipy^ 



Z X Tkj ~ , 
j^k I ^kj I 



(4) 



Since ip is single valued, the right hand side of the above equation has exactly the same multival- 
uedness as the left hand side. In other words we have, 



Uif. ^( 

i<j ^^'^^ 



a/2 



1 Zi , Zi 



n I 

i<j ^ ■' 



(5) 



^Later we replace a with Uq to denote the fact that this is the statistical parameter appearing in the quantum 
mechanics of anyons. While this is a dimensionless parameter here, it is not so in the classical case. This distinction 
is important in the discussion of semiclassics. 
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where 

Dkip = [Vfc + Ak]^p 

and 

z X fkj 



Akifk) = ^«y^ ■ 



This aUows us to write any higher order differential operators on ip in terms of corresponding 
covariant differential operators on the single-valued wave function -0. In particular a Hamiltonian 
operator, typically — V? + V can be written similarly. An eigenvalue equation written in 
terms of '0 can then be recast as a corresponding equation in terms of '0 involving the covariant 
derivative. 

Although both formulations are equivalent, dealing with operators on multivalued wave func- 
tions is much less transparent than dealing with operators on single- valued wave functions. Naive 
commutation rules, symmetries that one would expect by looking at an operator on single valued 
functions are not at all true in general for the "same" differential operators acting on multivalued 
wave function. 

Considering eigenvalue problem in terms of "0 has other advantages too. Since all the subtleties 
of multivaluedness are equivalently transcribed in terms of additional "interaction" terms (the so 
called statistical interactions), the eigenvalue problem is amenable to approximations. One is also 
on firmer ground in doing usual algebraic manipulations with operators. With these in mind we 
will work with single- valued wave functions with "statistical interactions" . 

As a first step one would like to understand the system of "free anyons". However, the 
statistical interaction falls off as | rjj as | rij |— > oo. So one is not sure whether the Hamiltonian 
with only statistical interactions has only discrete eigenvalues. One can put the system in a box 
to ensure discrete eigenvalues but then one needs suitable boundary condition. An oscillator 
potential ensures discrete spectrum without introducing a finite size. One could take some other 
confining potential but in the limit a ^ one should know the spectrum. One then has hope 
of doing at least the perturbative analysis [23] • Since the statistical interaction depends only on 
relative separations, the Centre of Mass (CM) dynamics should play a trivial role and oscillator 
potential also allows a separation of CM and relative coordinate dynamics. The oscillator potential 
problem can also be mapped on to a problem of anyons in a real, constant external magnetic field 
along the z-axis |24j. Bearing these facts in mind, we choose the oscillator potential without 
further justification. In order to derive the thermodynamic properties of a system of anyons there 
exist well defined methods of eliminating the dependence on the oscillator frequency |13| 120]. 

The Hamiltonian we consider is - after carrying out the usual scaling of variables - in terms 
of dimensionless quantities 



N N ^ p 2 ^ 

i=l i=l j>i=l *J i^jyk *J 

where 

= {fi -fj)x {pi -pj). 



and all distances have been expressed in units of y Notice that the statistical interaction is 
independent of the centre of mass coordinates. This is the Hamiltonian we analyze in the following 
sections both quantum mechanically and classically. In the quantum analysis, the Hamiltonian is 
considered to act on wave functions which vanish suitably at the coincident points. 
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3 Quantum Spectrum of Many Anyons 



In this section we summarize the known results on the quantum mechanics of N-anyons confined in 
an oscillator potential. A detailed discussion of the results may be found in the books by Lerda[3] 
and Kharejlj. The latter contains the results for several other types of confinement potentials as 
well. 

We first briefly summarize the results in the case of two anyons which is exactly solvable. In 
the case of N anyons {N > 3) the problem is not exactly solvable though there exists a subset 
of exact solutions. The results of this section provide the necessary framework to discuss further 
results in later sections. 

3.1 Spectrum of two-anyons 

The quantum Hamiltonian in the case of two anyons may be written as 

H = Hem + Hrel, (J) 

where Hem is the Hamiltonian that describes the dynamics of the centre-of-mass. Since the 
statistical interaction is translation invariant, the centre-of-mass part is independent of this. As 
a result the spectrum of this Hamiltonian is the same as the spectrum of a two dimensional 
oscillator. That is, 

Ecm = fuv[2ncm + |^cm| + !]• (8) 

Here ricm and Icm denote the radial and angular momentum quantum numbers of the centre-of- 
mass excitations. 

The relative Hamiltonian is given by, 

Hrei = f^[p'^ + -a-!^ + ^\], (9) 

where p is the momentum operator. The eigenvalue equation is easily solved by noting that the 
additional a dependent terms, may be combined with the centrifugal barrier with / shifted by 
I — a. The spectrum of the Hamiltonian is then given by, 

Erei = hw[2n + \l - a\ + I]. (10) 

Together with the spectrum of Hem, this provides the complete spectrum of two anyons in an 
oscillator potential. Indeed this is the only case that is solved fully. 

We next prove the existence of a subset of exact solutions when the number of anyons is more 
than two. Later we show that these solutions indeed correspond to the quantization of collective 
coordinates of the N-anyon system. We close the discussion with an analysis of the so called 
"missing states" by specifically taking three anyons as an example. 



3.2 Exact Solutions 

For discussing the known class of exact solutions it is convenient to use the complex coordinates 
Zi , Zi in terms of which the N-anyon Hamiltonian eq. © takes the form (here after we set h = c = 
1), 

H = -2Y^m + lY: -^--^ - « E(- - ^) + T E (11) 
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where = dij = di — dj, etc., and the eigenvalue equation is, 

H^{z„z,) = Eiljiz,,z,). (12) 

The conserved angular momentum J with eigenvalues denoted by j, is given by, 

J = Zidi - Zidi. (13) 

At large distances, the oscillator confinement dominates which immediately signals the pres- 
ence of the Gaussian factor in the wave function. Since the potential is singular as the relative 
distances go to zero, one needs to regulate the wave function in this limit to get the valid spectrum. 
We do this first by defining the variables: 

N N 

X = '[[zij; t = ^ZiZi, (14) 

i<j i=l 

We may classify the known exact solutions as follows: 

(a) j < solutions: 

= |Xr Mzi) e-*/2 ; j < (15) 
with the energy eigenvalues given by, 

E^N-i^a!^. (16) 

Here and in what follows 0j generically denotes an cigcnfunction of the angular momentum 
operator J with eigenvalue j. For elucidating the solutions the specific form of (f) is irrelevant. 

(b) j = solutions: 

V'o = \xr Mt) e"*/' ; 3 = (17) 

with i;^o(*)(= Ylk=i ^kt^) ^ polynomial of degree m in t. The corresponding energy eigenvalues 
are given by 

E = N + 2m + a ^^^~ . (18) 

The second solution is necessarily bosonic since t is symmetric, whereas the first solution needs 
explicit symmetrization and antisymmetrization of the wavefunction in terms of zi to obtain the 
bosonic and fermionic wavefunctions. Since this is always possible, the degeneracy of the first 
type of solution is exactly the same for both bosonic and fermionic type solutions for any given 
angular momentum j (< 0). We may also take a combination (product) of the solutions of the 
types discussed above, to get further j < solutions. This then generates an infinite tower of 
radial excitations for each value of j . 

(c) i > solutions: 

= </.,(zi) e-*/2 ; _^->o (19) 

with the energy eigenvalues given by, 

E = N^,-am^. (20) 

Caution must be exercised in choosing the value of j for these solutions since the wave function 
is not square- integr able for all values of j. In fact the lower bound is obtained by requiring that 
the wave function be square-integrable over the whole domain of a(0 < a < 1). This means that 
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j > {N — 1){N — 2)/2. If however this condition is not satisfied then the wave function remains 
regular only for some values of a (0 < a < 2j /N{N — \)) but not for all < a < 1 which gives rise 
to the so called non-interpolating solutions which have also been discussed in the liter ature|25j. 

All these solutions for the energy eigenvalues have a linear dependence on a with a coefficient 
zizN(N — l)/2 while the corresponding eigenfunctions are finite order polynomials apart from 
the overall | X and the Gaussian factors. These solutions (a)-(c) cover all the known exact 
solutions. However it is by now known that these exact solutions span only a subspace of the full 
Hilbert space and existence of nonlinear solutions has been shown numerically as well as through 
mean-field calculations. It is our aim here to understand the reason for the existence of this 
dichotomy. 



3.3 Missing states: Asymptotic analysis of the three-anyon problem 



We now discuss the question of "missing" states, a terminology used by Yong-Shi Wu|llj. The 
existence of these states has been shown numerically |12| for three and four anyons, and through 
perturbative analysis |131 114j for three anyon problem. Thomas-Fermi method has been used 
to obtain approximate ground state and excitations for large number of anyons |25- We focus 
on these states now. We cannot easily analyze these in general for a general N-anyon system. 
Therefore we restrict ourselves to the three anyon system, where it is well known that solutions 
which have nonlinear dependence on a exist. Even here, we can only comment on the nature 
of such solutions, but we do not see a simple method of obtaining even low lying excitations 
analytically. 



Recall that the two asymptotic regions, namely Zij 
Hence we may still write the full solution as. 



oo and Zij — > are always admissible. 



ij{zi,Zi) = exp 



^^ZiZij \Zij\°' tp2{Zi,Zi), 

i / i<j 



(21) 



with 1^2 satisfying the equation. 



-2 m + YizA + m) - 2a ^ 



1p2 



E-N 



a 



N{N -I) 



■02 



(22) 



The known exact solutions for the energy eigenvalues, discussed previously, have a linear de- 
pendence on the statistical parameter a while the corresponding eigenfunctions are finite order 
polynomials, ^/'2) apart from the | X |" and the Gaussian factor. In fact a simple scaling argument 
shows that if is a polynomial (i.e. has a finite degree), then the corresponding eigenvalue must 
be linear in a. For, if -02 is a polynomial in Zj, Zj, with highest total degree d then for Zj, — > oo 
the polynomial becomes a monomial and only the scaling operator term, dominates. This 
gives, 

N{N -I) 



E = N + d + a- 



(23) 



Thus to have non-linear dependence on a, ■02 can not be a polynomial. 

Suppose 02 admits a power series representation with infinite radius of convergence but the 
series does not truncate then the scaling argument fails. However the analysis of such a series 
solution always seems to lead to exponentially divergent behaviour making the solution non- 
normalizable, ie., the power series has to truncate. But then only linear solutions are possible 
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because of the scaling argument above. The other possibiUty then is that the power series has a 
finite radius of convergence, i.e., a "scale" R has to enter if nonlinearity is to be possible. 

Let us pursue this scaling argument further, consider ■i/;2(Azj, Azj) where A is the scaling pa- 
rameter. Having a scale R, radius of convergence, means that ip2 has two different representations 
as the scale parameter A ^ and A — oo, each being valid for A < Xmax and A > Xmin respec- 
tively. For A > Xmin the scaling argument can still work but now the series need not have only 
integer powers, i.e., d can be a nonlinear function of a. We conclude then that if a nonlinearly 
interpolating eigenvalue is possible at all, the corresponding eigenfunction must have two different 
series representations as A — s- and A oo. One then has to try to match the two series suitably. 

Exploring this possibility is extremely complicated for general N. The first non-trivial case 
of = 3 is analyzed in [21] . While this analysis shows that the possibility of finite radius of 
convergence is viable, one is not able to obtain any specific non-linearly interpolating solution. 
Observe that to get a discrete spectrum, the solutions of the eigenvalue equation must be regular 
and normalizable. For a single power series, this forces us to truncate the series and one obtains 
linearly interpolating spectrum. For multiple power series for the eigenfunctions, while non-linear 
dependence is conceivable, it is hard to see regularity /normalizability to get discrete spectrum. 
One can get a glimpse of the intricacies through the following analysis. 

Consider the method of analyzing solutions in the three anyon problem consists of using the 
Fourier expansion of the wavefunction. Since the center of mass degree of freedom is irrelevant 
for further arguments, we may separate that out and write the wavefunction as 



i/j = exp 



where ri,r2 are relative coordinates chosen appropriately. In general (j) may be written as, 

Mn,r2,eue2) = Yl Xn„n.(ri,r2) = ^e^^'^e^^ X,,n(ri,r2) (24) 



ni,n2 



for a given angular momentum j, where j = ni + n2, n = rii — n2. We can now use this 
representation to look at the solutions for a given j. The eigenvalue equation now becomes 
Hrei4'j = ^'Pji where ff^ei is the three-anyon Hamiltonian in relative coordinates and £ = E—2—3a 
as before. Substituting for (pj the eigenvalue equation takes the form, 



3xo 



92 

42;i--^ + 4{l + a- xi)—- 
oxf 



d_ 

' dxi 



{j + n)[3 + n - 4a) 



+ 



q2 q 

4x2— + 4(l-h2a-X2) 



4x1 

(j - n- 8a) 



' dxf, 



Xl 



4:Xi—^ + A{l + 3a - Xl 
oxf 



+ 



X2] 



_d_ 

dx-2 



dx2 

d 

dxi 
{j - n 



4x0 



+ 2£ 



Xj,r. 



{j + n + 4)(j + n + A- 12a) 



if 



4X2 



+ 2£ 



4x1 



where x,- 



This expression is exact, representing an infinite set of coupled equations. They 



relate Fourier modes differing by 4, i.e., if xo, Xi,X2,X3 are known then XAk,X4k+i,X4k+2, XAk+3 get 
determined in terms of xo; •••) X3 etc. However there is no relation among xo, X1tX2, X3 themselves. 
In a sense these four functions will give four independent solutions of the eigenvalue equations. 
We can then deal with a given "tower" separately and independently and this is true for every 
given j. Let us concentrate on one tower. Now three distinct cases arise naturally: 
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(a) Only one member of the tower is nonzero, ie., Xj,n = Xj,m^n,m- 

(b) Only a finite number of Xj,n^ a^re nonzero, ie., Xj,n = V n > ni and V n < n2 with 
n2 < rii. 

(c) Infinitely many Xj.n's are nonzero. This gives raise to the following cases: 

Xj,n = V n > ni 
Xj,n = V n < n2 
Xj,n / V n. 

Case (a) is simple to analyze and it reproduces the known exact solutions. In the exact 
eigenvalue equation relating Xj,n 

and Xj,n+A given above, we set the left hand side to zero which 
gives an equation for the Xj,n- But the right hand side also gives rise to another equation for the 
same Xj,n when the coupling to Xj,n-i is taken into account. However the form of the equations 
suggests that the two equations are separable but must be consistently solved. The consistency 
conditions immediately yield the linear exact solutions which are already outlined above. Cases 
(b) and (c) seem quite complicated and nonlinearly interpolating states must be in one of these 
cases. Consider the normalization condition on the full wavefunction, 

II f= E / dx^dx2\X\^''\xj,n?e--'+-' = y^Cn. (25) 

n 

If V'j is given by case (b) then there are only a finite number of terms in the norm and if CnS, 
are finite then '^j is normalizable. Thus the quantization condition for the nonlinear states must 
arise by the demand that Cn must be finite. For the case (c) even if all C^'s are finite we may 
still get the norm of the wavefunction to be infinite because the sum above may not converge and 
then the quantization condition would be the convergence of C„. It is therefore apparent that 
non-linear dependence on a must come about in a very complicated manner. 



4 Classical Analysis of Many Anyons and Integrability 

Having noted the difficulties inherent in describing the non-linearly interpolating solutions, it is 
striking that the subspace spanned by the exact solutions are obtained so easily. Obviously these 
two sets must have very different origins. Indeed as we will see below, the exact solutions have 
their origin in the quantization of collective degrees of freedom. To explore this aspect we analyze 

the system classically. 

We may now, therefore, ask the question, what is the classical Lagrangian which gives raise 
to the Hamiltonian of Eq.(6)? In fact it is quite straight forward to see that the Lagrangian is 
given by, 

1 ^ ^ / _ \ 

^ = 9 E ^^^f + " E - ^(^"'^); = !^ , (26) 

^ i=l i<j ^3 J 

where denote particle coordinates and V(fi) is some confining potential which we choose to 
be the harmonic oscillator potential V(fi) = ^^i^^j^- The a-dependent term is the statistical 
interaction. The harmonic oscillator potential is convenient since the dynamics of the system is 
well understood in the limit a = 0. As mentioned in the previous section, the corresponding 
quantum system displays the existence of a subset of exact solutions, though the full spectrum 
is not analytically solvable for many anyons. It is possible that the system is at least partially 
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separable to allow for the known exact solutions. We will return to this question to the next sub- 
section and consider the question of integrability first. If one considers the classical Lagrangian, 
the statistical interaction is a total time derivative. The Euler-Lagrange equations of motion are 
therefore the same as those of a 2N-dimensional oscillator and this of course is an integrable system. 
However the quantum spectrum displays not only the exact solutions but it has been demonstrated 
numerically in the case of three and four anyons that the spectrum also shows evidence of many 
level crossings and level repulsions jl21 I18j which is often taken to be an indication of non- 
integrability. Here we argue that even though we have 2N constants of motion in involution the 
invariant surfaces do not have the topology of a 2N dimensional torus. 

We begin with the analysis of the classical Lagrangian given by Ea.()26|). It is convenient to 
write the Lagrangian in the form, 



L 



1 ^ 



N 



(27) 



where the dots indicate the time derivatives. The first step is to introduce relative coordinates, 
also called Jacobi coordinates and separate the trivial centre of mass degree of freedom. To this 
end we write, 



Pa 



with the inverse relation given by. 



a = l,...,N - 1; 



(28) 



-. N-l 

I — \ ^ \ - 

^—Pi-l+ > 

7 ^ ^ 



Pk 



+ Rem = Pa + Rc 



(29) 



where pa, a = 1, N — 1 are dimensionless relative coordinates and Rem is the centre of mass 
coordinate. It follows that. 



It is straight forward to see that, 
where 



Lirel — 2 



L = LcM + Lrel, 
LCM = -j[RcM - R-Cm]^ 



^2 -Q 
Pa -Pa 



A\. Pa X 'p^ 
A'r, At p,.pd ' 



where Af- = A^ — A^-. From now on we concentrate only on the L^ei and drop the subscript. It 
is easy to see that the Euler-Lagrange equations of motion are. 

Pa = -Pa- (30) 



There is no a dependence in these equations since the corresponding term is a total derivative. 
From these equations it follows that 



1 



^a = l^[Pa + Pa\^ L = Pa X Pa] a = l,...,N -1 



(31) 
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are constants of motion. In the Hamiltonian formulation the pa are expressed in terms of conjugate 
momenta and coordinates which do contain the a dependence, 

dL _ ^ -^ij ^ij 

dpa. f^. Ai^ Af^ p,.pa 

P - _ ■ , 4 /I _ ^fj A\j pbx 

Opay A.J A.J p^.pd 

and the constants of motion in relative coordinates are, 

£a = ]^[{Pax + aAaxf + {Pay - aAayf + Pa"^] 

la = Pa^ Pa - a{PaxAay + PayAax)- (34) 

Clearly the Hamiltonian is given by 

N-l 
a=l 

Note that the o-dependent term in L is well defined only if rij ^ for all i, j, that is, only if 

A%A% pa.pb + V 

Consequently the expressions for fa, /a and H are also valid only if rjj is not zero. In effect the 
classical configuration space on which L is well defined is the space 

gAf-i ^ _ ^. A = { Pa / A\jA% p,.pd = for some pair(s) (35) 

The space Q^~^ is not simply connected. Its fundamental group vri is the same as the fundamental 
group of {R2 — {N — 1)} points which is known to be nontrivial. For N = 2, iix = Z while for 
N >2>, TTi is non-Abelian. 

The corresponding phase space is the cotangent bundle of Q on which £a-:la and H are well 
defined. This phase space is topologically nontrivial and for our purposes we do not need the full 
machinery for handling this topologically nontrivial phase space. It is sufficient to pretend that 
the full space is still R2 ^^ but simply avoid coincident points. 

It is straightforward then to prove the following Poisson bracket relations, 

{£a,£b} = {laM = {£aM = ^ V a, 6 = 1, iV - 1, (36) 
{H,£a} = {F,/4 = 0. (37) 

Thus for 4(A^ — 1) degrees of freedom P^, pa we have 2(A^ — 1) constants of motion in involution. 
So a necessary condition for a system to be integrable is satisfied. We will tentatively refer to 
this system as being integrable (or potentially integrable) in the "Liouville sense" . However, for 
the system to be integrable via action angle variables, further conditions have to be satisfied (see 
page 271 of [17]): If M(£a,la) denote the set of points in phase space at which the constants 
of motion have values £a,la then (a) M is a 2(A^ — 1) dimensional submanifold iff £a,la are all 
independent, (b) if M is a submanifold which is compact and connected then M is a 2(A^ — 1) 
dimensional torus. If these conditions are satisfied then one can introduce action-angle variables 
in the neighbourhood of M. The orbits on M will all be conditionally periodic in general and for 
the periodic orbits one may use Bohr-Sommerfeld quantization to get a subset of eigenvalues. 
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We will show in the next section that there exist M{Ea,la) for which action-angle variables 
can be introduced. These correspond to the collective motion of N-anyons and their semiclassical 
quantization reproduces the exact eigenvalues of the quantum system that are already known. 
However, if we go away from these initial conditions some of the orbits of the oscillator will cease 
to be periodic when a dependence is introduced. For instance the Euler-Lagrange equations 



For any pair ij, i ^ j, fij in general describes an ellipse in the configuration space. However for 
zero orbital angular momentum, the ellipse degenerates to a straight line, hence the point fij = 
is on such an orbit. Since these points are not admissible when a ^ these orbits cannot lead to 
periodic orbits in the phase space. For all choices of £a,^a such that for some initial conditions 
the phase space orbits will have rij approaching arbitrarily close to zero, M{£a,la) cannot be a 
torus. (Even if M is a manifold it will fail to be compact and/or connected.) In appendix A, 
we consider the explicit example of an angular momentum zero orbit in the case of two anyons, 
which necessarily passes through the diagonal point. Even though in configuration space this 
may appear to be periodic, in the phase-space the orbits stay on the surface of a cylinder which 
can be broken into two disconnected sheets. Therefore for a 7^ 0, although we have a potentially 
integrable system it is not generically integrable via action-angle variables. The existence of such 
an M indicates the possibility of "extreme sensitivity" to the initial conditions. This is a possible 
reason for the "level repulsions" seen numerically. This bears a resemblance to the billiard system 
considered by Richens and Berry 17 though the reasons for the failure of integrability via action- 
angle coordinates are different. Following Richens and Berry we conclude that our system is 
Pseudo-integrable. We will sharpen and formalize this notion later in Section 5. 

4.1 Collective modes and relation to exact states 

We may now ask the question, what causes these states, though a subset of the full Hilbert space, 
to be exact solutions of the N-anyon Hamiltonian which is in general non-separable. Indeed as 
we will show now, these solutions may be obtained by quantizing the collective coordinates of 
the N-anyon system in relative coordinates consistent with the equations of motions for a set of 
initial conditions. In some sense these collective coordinates do not give us any information on the 
internal dynamics of the system. These are therefore trivial solutions dynamically. The anyonic 
interaction term in L is invariant under two sets of time independent transformations: 



where U{9) denotes a rotation of the vector by an angle 6 in two dimensions. Therefore for any 
given configuration {pa} at any given time t we can always rotate the axes so that piy = (say). 
Thus by making a time dependent rotation we can ensure that rjiy = for all t. Defining, 



fi = — fi =^ fij = — fij (independent of a). 



Va = U{6)pa V a, 
Va = >^Pa V a. 



(38) 
(39) 




(40) 



The following identities are easy to prove: 



Pa-Pb = Va-Vb, 
Pa'X Pb = fJaX fjb, 

PaX Pb = 9f]a.f]b + f]aX fff,. 



(41) 

(42) 
(43) 
(44) 
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Now the set {ija; a — 1, A'" — 1} is effectively a 2{N — 1) — 1 dimensional vector since rji^y — Ofor 
all time t. We can therefore introduce the standard spherical coordinates by the usual procedure, 



0; 



1, 



where 



^N-a,x = SiS2...S2a-2C2a-l, 
^N-a,y = SiS2...S2a-lC2a, 
^l,x = SIS2...S2N-4, 



(45) 



(46) 
(47) 
(48) 
(49) 



where o = 2, A'^ — 1 and S/^ = sin^^, = cos^^. In terms of these variables the Lagrangian 
may be rewritten as 



L 



1 
2 

+ a 



i<j 



(50) 



Since fja are obtained from pa by the same rotation matrix for all a, the angle 9{t) clearly 
describes a collective rotation of all the N anyons about some reference axis. The anyonic interac- 
tion term ( a total time derivative) is manifestly independent of i?,(t). This may also be regarded 
as a collective mode as discussed below. It is the semiclassical quantization of these two modes 
that yields the exactly known energy eigenvalues. To elaborate these points let us consider the 
Euler-Lagrange equations of motion. Since the a dependent part of the Lagrangian is a total 
time derivative we can ignore it for analyzing the classical equations of motion. Clearly these 
equations are identical to that of the oscillator equations of motions. Translating this into the 
spherical coordinates we obtain, 



[R + R- Re^j^a - [RO + 2R6]VU + 2[i? - ReV]ia + R^a = 0, 



(51) 



where 



V = 



1 
-1 



and in the matrix equation above the ^ is a column vector with two elements. The matrix V 
between two vectors essentially generates their cross product. Taking the dot product with 
and summing over a we find. 



[R + R- RO^] = 2Re ^ ia X ia + RY.^a, 

a a 

where we have made use of the identities, 

^ ] ^a-^a — 0) ^ ] ^a-^a — ~ ^ ] ^a- 



(52) 
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Taking the cross product with and summing over a we find, 

[Re + 2^^] = -2R ia xia-RYl ia X Ca- (53) 

a a 

Using the above two relations Ea. H51|) can be rewritten as, 

RCa + 2[R - ReV]ia + [2R9ib x ib + R^b]ia + mib X ib + RCb X ibWCa = 0, (54) 

where sum over b is assumed. Here the equations of motion for R and 9 are stih coupled to all 
the other internal coordinates, and as yet they are not collective coordinates. We therefore 
need to impose a set of initial conditions which may lead to the separation of these two modes 
as collective modes from the internal variables 0^. To realize this let at time t = 0, all velocities 

^a(t = 0) = for all a. Then the above equations reduce to, 

R[ia + (5Z ib X ^bWia] = 0, (55) 
b 

R + R- Re"^ = 0, (56) 

RO + 2Re = -i? ^ i; X ^b- (57) 

b 

Now consider the first equation for i? ^ and a = 1. From the initial conditions it is obvious 
that, 

iix = 0; 6x J2^ib X ib) = 

b 

because S.iy = V t. 

Now if 6.1X is nonzero, then Xlfe'^b x ^b = 0, and hence 

6 = V b. 

Since the equations of motion are second order in time t, we have proved that: if 6(0) = 
V a, and R{0), 6x(0) are non zero thenM t, 

6(t) = 6(0), 

R + R- Re"^ = 0, 
Re + 2Re = 0. 

For oscillator (a = 0) R(0) and 6x(0) being non zero is a special class of initial condition. Indeed 
in general for a second order equation system if both first and second order derivatives vanish at 
any t then this can hold only for some restricted class of coordinate values. For the equations of 
motions considered in the configuration space R2~^ the restricted class is precisely characterized 
by i?(0) 7^ and 6x(0) / 0. However, when a ^ the configurations space is Q^~^ which does 
not contain points which have R = and 6a; = 0. Thus there are no restrictions on the initial 
conditions in Q^~^. 

The above initial conditions amount to freezing the "internal motion" of the anyons. The 
remaining motion is a collective motion described by R(t) and e{t) which is just a 2(A^ — 1) 
dimensional oscillator. R(t) being nonzero implies that the angular momentum must be nonzero. 
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This explains in what sense R(t) can be interpreted as a collective mode. For describing the 
motion of collective mode the effective Lagrangian is 

Le// = \[R' + RH^ - R'] + a^^^^^e. (58) 

This is identical to the Lagrangian in the relative coordinate for a two anyon system with a — > 
^ Af(JV-i) 2('jY _ dimensions. Semiclassical quantization will then reproduce all exactly 

known energy eigenvalues summarized in Sect. 3. In this effective Lagrangian picture the only 
known memory of N resides in the coefficient of a. This can be understood by noting that when 
all the N-anyons are rotated by 27r about the centre of mass they also circle each other to pick up 
the extra phase. 

If on the other hand we consider i? = ^ = Oatt = class of initial conditions then we see that 
R,9 arc nonzero and hence R{t),6(t) do depend on 0^. Thus the collective motion is not fully 
decoupled from the "internal motion" . We therefore refer to this system as partially separable. 

We also note that the same conclusions may be drawn from the Hamiltonian formulation. For 
completeness we give the relevant expressions and arguments. The conjugate momenta for the 
hyper-spherical variables are given by, 

Pr = R, (59) 
Pe = R'e + R'Y.^,F, + a ^^^- ^\ (60) 

P^ = R\le^ + R^eF^ + aG^, (61) 

where, 

Tlx = siS2...Si^-i, tJ, = l,...,2N -3, ri = 1 
and Ffj^ and Gfj, are defined through, 

2Ar-4 JV-1 

E ^^^^ ^ E ^« X £ 

/i=l a=l 
fM=l i<j ^ij ^ij ?c-?d 

The Hamiltonian is then given by 



where 



H = Hi + H2, (62) 



= + + (63) 



H2 = 



, N(N-1) ^ F^(P^-aG^) ^2 



1 

2B^ 



Ep2 

1^ ^ 



(64) 



It is easy to see that for the special initial conditions 9^ = 0, H2 = 0. Also the Poisson bracket 
{Hi, H2} (X H2 vanishes for these initial conditions. The corresponding quantum statement would 
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[Hi, H2] oc H2, and therefore if we consider the subspace of states {V^} on which H21P = then Hi 
will act invariantly on this subspace. Further H = Hi on this subspace and thus the eigenstates 
of Hi will be exact eigenstates of the full system and conversely. But the problem of solving Hi is 
analogous to that of two anyon problem with a —>■ a ^*-^ Therefore these solutions are given 
by, 

E = 2m+\j-a ^^^~'^\ + {N-l), (65) 

4, = C7i?IJ-"le-^'/V(^), (66) 
Hiijj = E^P (67) 

where f{R) is some polynomial of degree 2m in R. In general the normalization constant C 
in ^ will have dependence on 0^'s, restricted by H21P = 0. These are precisely all the known 
solutions and the above argument shows that there are no more solutions satisfying Hiip = Eip 
and H21P = 0. Thus H21P = characterizes the subspace spanned by all the known exact solutions. 
So the quantum counterpart of classical partial separability implies the existence of these solutions. 



5 Semiclassical Analysis — Periodic Orbit Theory 

We have seen the difficulty in obtaining the non-linearly interpolating eigenvalues. We have also 
discussed their possible relation to, at classical level, pseudo-integrability on the one hand and the 
possibility of two different asymptotic series at the quantum level on the other hand. There is yet 
another possible approach to gain some understanding of these eigenvalues, namely semi-classical 
analysis. This is what we discuss now. We begin by classifying the classical orbits for N-anyon 
system and discuss the trace formula in the subsequent section. Most of the discussion of this 
section is based on reference |28j . 



5.1 Classification of orbits and symmetry reduction 

We will first carry out a classification of orbits for anyons and see how the pseudo-integrability 
amounts to a reduction of the dynamical symmetry. We will also see the role played by the 
fundamental group of the phase space in this regard. Noting that a trajectory is specified by 
giving an initial point, we will use the dynamical symmetry of the oscillator system to group 
various orbits into continuous families. For anyons, certain points and hence certain orbits are 
disallowed (removed) with the result that one gets several continuous families. 



5.1.1 The case of two anyons 

As usual we can separate the center-of-mass (CM) dynamics from the relative coordinate dynam- 
ics. The CM dynamics is trivial - identical to an oscillator locally and globally. The anyonic 
feature is contained entirely in the relative dynamics described as: 

Q = 

L = i(r^-r^)+ac^ 

H = ^{TT^ + f^) 

The classical trajectories are trivially known in these cases. These are ellipses with two 
constants of motion, the energy (E) and the angular momentum (J = (r x r)z)- The sign of 
J gives the sense of traversal and a trajectory reaches the point r = if and only if J = 0. 



r = Q X i?2 

e = tan-^y/x (68) 
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These degenerate trajectories ( J = 0) are not allowed for the anyons though are allowed for the 

oscillator. 

Consider the set of all possible trajectories with a given energy. Fix two such trajectories with 
angular momenta J, J'. Suppose there exists a continuous family of interpolating trajectories 
connecting these two. Clearly if the angular momenta have opposite signs, then at least one of 
the interpolating trajectories must be degenerate. But this is disallowed for anyons. Hence, for 
anyons there must be at least two 'orientation' classes of trajectories distinguished only by the 
sign of angular momenta. 

While we have shown that there must be at least two families we do not know yet if there 
must be precisely two families. We will use the symmetries to show that there are precisely two 
families for anyons and precisely one family for the oscillator. 

As far as the PB algebra is concerned, the oscillator and anyons are identical and hence we 
do have infinitesimal symmetries forming the Lie algebra of ^7(2) (or 0Sp{4,R)) for both cases. 
The constant H surface {S^ in the phase space for oscillator) is of course invariant under the 
infinitesimal symmetry transformations. For anyons, the phase space has r = removed and the 
constant H surface has a great circle (one dimensional) removed from the S^, say a "longitude". 

In the appendix B we have given the details of the OSp{4N, R) symmetry of the oscillator. 
In the present case of two anyons with center of mass coordinate separated, we have effectively, 
A'^ = 1 and the symmetry generators form OSpiA^R). We have only the 4 generators of the Tj 
type. Of these four, ui generates rotations of position and momentum while U2 generates time 
evolution. The generic integral curves are given in the appendix B. 

From appendix B we also recall that 05*^(4, R) acts transitively on H = E sphere. This shows 
that for the oscillator there is precisely ONE family of trajectories. For anyons we need analogous 
result with the extra condition that degenerate trajectories are avoided to conclude that there are 
precisely two families. 

Claim : Given any two trajectories with the same sign of their angular momenta, there exists 
a one parameter subgroup connecting the two without changing the sign of angular momenta. 

Proof : Observe that Ui,U2 transformations leave the angular momentum invariant. Choose 
u = cos(/?)u3 + sin(/?)u4. The one parameter group generated by this generator transforms 
the angular momentum as, (appendix B) 

2J{a) = 2J-sin^(cr){2J-a;uLua;} + sin(2cr){a)Lua;} (69) 

For the choice of u made, uLu = — L and putting = 

J{u) = J cos(2o-) + sin(2c7) 

= ^J'^ + Jl cos(2a-5) , b = i&a-^{Ju/J) G (-7r/2, 7r/2). (70) 

This is valid for all choices of (3. Observe that Ju depends on (3. Furthermore it satisfies. 



3(3' 



2 



r) T 

MP) = Jn(0)cos(/3) + ^(0)sin(/5) ^ (71) 

{Ju)max = { {Ju{0)f + (^(0))' y = V - J'' 
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The last equality follows by explicit evaluation. Choosing /? to maximize Ju then implies that 
the square root in the equation for J(ct) is just the energy E. Since for any given energy, we must 
have J'^ < E'^, it follows that J{cr) covers the entire range of possible J values. 

From this it follows that we can always choose a a such that J(ct) = J' without passing 
through zero angular momentum. This proves the claim. Hence there are precisely two families of 
trajectories each spanning the whole orientation class. This gives the classification of trajectories. 

One may view the S"^ as compactification of with the south pole as the origin of i?^ and the 
north pole as the point at infinity. The removal of A now corresponds to removal of an infinite line, 
say, the z-axis of the R^. It is immediate that the fundamental group is the group of integers, Z, 
and therefore the basic trajectories winding around the z-axis belong to two orientation classes. In 
this case the fundamental group directly shows how the basic trajectories are naturally classified 
into the "orientation classes". 

Consider now the Hamiltonian vector fields generating the three one-parameter subgroups of 
SU {2), ui generates rotations about the z-axis. Its orbit through any generic point avoids the 
z-axis and the vector field remains complete (its integral curves range over the full real line). 
The other two subgroups on the other hand have orbits necessarily cutting through the z-axis. 
These vector fields therefore are necessarily incomplete (for the anyon case). The transformations 
generated by these vector fields do not form a group and the dynamical symmetry group for 
anyons is reduced from U{2) to U{1) x U{1). 

Thus the removal of the set A has altered the fundamental group, has classified trajectories into 
orientation classes and has also reduced the dynamical symmetry group form U{2) to U{1) x U{1) 
due to the incompleteness of vector fields. 

Note that while the symmetry group is reduced from U{2) to U{1) x ^7(1) for anyons, the 
rank has remained the same and therefore the two-anyon system continues to be integrable. The 
integral curves being circles, we have integrability via action-angle variables and we get the exact 
spectrum from the semiclassical approximation. In fact one can follows the EBK quantization 
procedure to reproduce the exact spectrum. 

5.1.2 The case of many anyons 

The N > 2 case differs from the N = 2 case in many ways. The potential symmetry group, 
OSp{4:N, R), (see appendix B) generating families of trajectories is lot more complicated and so 
is the fundamental group for the phase space of anyons. However the analysis of the previous 
section already suggests a strategy for obtaining a classification of the families of trajectories. 
A generic trajectory may be viewed as a collection of ellipses traced by the fij,i < j. 

To get a convenient handle on the disallowed trajectories, define Jij = {fij x fij)z- Clearly a 
trajectory will cut through the set A if and only if vanishes for at least one pair of indices. 
We will term such a trajectory as degenerate. For anyons only non-degenerate trajectories are 
allowed. Note that Jij are constants of motion just as Jj and Ei = i(r? + pf) ||16j are. Each 
of these elliptical trajectories will have a sense of traversal ( Cij = sign of Jij ). Thus a basic 
trajectory of particles has an associated set: {fiij}. Exactly as in the previous sub-section, 
the set of non-degenerate trajectories can be classified by these 'orientations'. Provided we can 
exhibit trajectories which realize all possible choices of eij, we will have at least 2^(^~^)/2 families 
for anyons and precisely as many if every member of a class is connected to every other member 
of the same class by a symmetry transformation without crossing A. It is easy to show that every 
choice of Cij is realized by considering concentric circular trajectories for each of the particles 
with various possible signs for the angular momenta Jj and various possible ordering of the radii 
of the circular trajectories. For oscillator we will have precisely one family if OSp{4N, R) acts 
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transitively on the constant energy sphere, S^^~^. 

The proof of transitivity in the case of oscillator is most directly obtained by noting that 
OSp{AN, R) is isomorphic to U {2N, C) and it is well known (and easy to see) that the coset 
space U{2N,C)/U{2N - 1,C) is difFeomorphic to S^^-'^. This immediately shows that for the 
oscillator, there is precisely one family of periodic trajectories. The counting of families is more 
tricky in the case of the anyons. Let us see the symmetry reduction in a slightly different manner. 

Let A' denote the set of all degenerate trajectories i.e. at least one of the Jjj's is zero. Note 
that A and A' are different, the former being a proper subset of the latter. Suppose we find the 
subgroup of symmetries which leave A' invariant, then the same subgroup will also leave the set 
of all non-degenerate trajectories invariant. Since we are looking for a group, we can consider 
the infinitesimal action. Let J = Yl^^jJij- Clearly, J = characterizes the set A'. If we 
have a degenerate trajectory with two or more Jjj 's zero, then every infinitesimal action will keep 
within the set of degenerate trajectories. If only one Jmn is zero then infinitesimal action will 
act invariantly if and only if J„j„ is invariant. Of course we could have degenerate trajectories 
with different angular momenta being zero and hence for the invariance of A' it is necessary 
that the generators of the subgroup must leave each of the Jij 's invariant whenever Jij itself 
is zero. This subgroup is identified in the appendix B. There are two forms of generators. Any 
generator with all the block matrices being the same u matrix with u G OSp{A, R) will leave all 
the J's invariant. These transformations act on the center-of-mass variable alone and constitute 
the expected OS'p(4, R) symmetry. In addition, generators with only diagonal blocks being the 
same u and all the off-diagonal blocks being also leave the J's invariant provided u = Ui or 
U2 (or a linear combination thereof). These are just the total Hamiltonian and the total angular 
momentum - Ji) which generate time evolutions and common rotations of all the positions 
and momenta. These are the only subgroups leaving the set of degenerate (and non-degenerate) 
trajectories invariant. The vector fields corresponding only to these will remain complete. Thus 
wc sec that the symmetry is reduced to OSp{A, R) x 0(2) x 0(2). This is identical to the symmetry 
in the case of two anyons. 

It is apparent from the above discussion that pseudo-integrability of anyons really means 
that not all the Lie algebra level symmetries exponentiate to Lie group symmetries. One should 
distinguish between "infinitesimal integrability" and integrability. A general discussion is given 
in the sub-section 5.3 . 

This purely classical analysis has already explained the qualitative result that for anyons only 
the total energy and the total angular momentum are the conserved quantities (good quantum 

numbers). Since wc also have a classification of periodic trajectories wc will take the next step of 
attempting semiclassical approximation to see what further information can be obtained. 

5.2 Semiclassical Spectrum of anyons 

To develop a semiclassical approximation the central quantity of interest is the propagator defined 
as, 

= J](E-E„ + !e)-' (72) 

n 

= ^ riiE-E^)-')-i7rY, KE-En). 

n n 

In the last equality, V denotes the principle value while the second term is of course the density 
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of states. Defining G = iG allows us to write the density of states, g{E) as , 

g{E) = l-{G{E + ie) + G*{E-ie)} 

= - Re{G{E + ie) } (73) 
vr 

The derivation of the semiclassical trace formula begins with the definition of the propagator, 
G{E), with trace operation expressed via a path integral representation. The trace operation 
combined with the stationary phase approximation (SPA) gives the propagator as a sum over 
periodic orbits (in the phase space), of terms with an amplitude given by a suitable Van-Vleck 
determinant times a phase whose exponent is the value of the 'action', J p • dq around the periodic 
orbit together with a contribution from the Maslov indices. This works very well for systems 
having only isolated periodic orbits. For systems with symmetries, as for the oscillator or anyons, 
the periodic orbits come in continuous families and a modification is needed. 

In cases with symmetries, the sum over orbits gets replaced by a sum over families of orbits 
together with a measure factor. If the family arises due to a symmetry, as is usually the case, 
the amplitude and the phase is same for all members of the family and as such can be computed 
from any one member. This is discussed by Littlejohn et al [20] in detail. The origin of families 
of periodic orbits in the case of oscillator is of course the 05p(4, R) dynamical symmetry and 
for anyons essentially the same transformations generate families of trajectories. This is already 
discussed in the previous sub-section. 

A system could also admit the so called "diffr active" orbits. Generically these are orbits which 
are closed in the configuration space but not in the phase space (or that the Hamiltonian flow is 
discontinuous) 30 . Within the framework of so called "uniformly continuous approximation" |31j . 
these lead to further contributions to the semiclassical propagator. The anyon system presents us 
with both these features . 

In the context of semiclassical approximation, one point needs to be noted. Observe that in 
terms of the tTj , fj coordinates there is no explicit reference to a and A is also precisely the set 
on which one gets (5- function singular Poisson brackets with a appearing as the coefficient. If 
one did inverse Legendre transformation to get a Lagrangian one would not get the a dependent 
term, but of course the configuration space will have the A removed. Using pijfi variables and 
then doing inverse Legendre transformation will produce the Lagrangian with the a dependent 
total derivative term. Now the removal of A is made explicit by the a dependent term. It is this 
explicit form which is convenient for computation of action for periodic orbits. 

The action for the oscillator is equal to {27rE)/uj and the action for anyons is the same as that 
for the oscillator except for the contribution from the total derivative term in the Lagrangian. 
Its value is ib27ra with sign depending on the orientation class(es). We have already seen that 
for = 2 we have precisely 2 orientation classes while for general A^ we have at least 2^(^-i)/2 
classes. 

Are there "diffractive contributions"? To appreciate this, let us consider the two anyon case 
in more detail. The propagator can be computed exactly from the known exact spectrum |32j . 
The exact spectrum is given by, 

En,j = nuj{2n + \j -ag\ + l), n > 0, j e Z. (74) 
From this the partition function is obtained as, 

7fr3\ _ cosh{ Phw {ag - 1) ) + cosh( f3fuvaq ) 

^^^^ - 2sinh\l3fu.) ' ^^^^ 
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and the exact density of states is obtained as: 



g{E) ^ — / e^^ZiP) 



E , v-^, .27rkE , , ,2iTkE , 

|1 + / ,iCos(— h 2'Kkaq) + cos(— 2TTkaq)\ 



^ ' k>l 

1 r / \ ,2TTkE . , 

Z^i^Og sm(27rA;ag) sm(-^)| (76) 



fe>i 

^ ' k>i 

Rewriting the products of sines as differences of cosines, we get 

E 



9{E) 



1 , E .2TTkE , , 

+ ?^( cos(^ + 2vrto,) 

k>l 

I , E ,2'KkE , , 
+ ^ ( ^ - "g ) 27rA;a5) (77) 

k>l 



L T, (-1) 



k>l 



cos(^: 1- Trkuq 



2^ (-1) cos(^ vr/ca 

A:>1 



Remarks: 

1. The partition function is manifestly invariant under — > 1 — Og . The convergence of 
integrals in computing the density of states for Og £ (0, 1) requires ^ to be greater than or equal 
to 1. 

2. The first term above is the usual Thomas-Fermi term [HSj. This will be suppressed in 
expressions below. 

3. In the limit ^ — > keeping E^ Uq fixed, the exact propagator effectively looses all dependence 
on aq. This is also the semiclassical trace formula for the oscillator showing that the semiclassical 
approximation is exact for the oscillator. All the other terms are sub-leading. 

4. For a comparison with semiclassical approximation, a different limit is implied. One should 
use Oc = haq and the SPA will be done keeping Oc fixed. Then only the last two terms will be 
sub-leading and the propagator will continue to have Uc dependence. This is the limit considered 
in the following. 

5. In the last two sub- leading terms the argument of cosine is half of that in second and the 
third term. As is well known, the argument of cosines is the classical action, J p ■ df, around 
a classical orbit including possible contributions from "Maslov indices". The last two terms are 
therefore suggestive of contribution from "half orbits" . Indeed if one takes a Fourier transform 
of g{E) ^3 , then one sees a peak at 1/2 period from the last two terms. One may therefore 
suspect a contribution from classical trajectories with half the period. Do we have such classical 
trajectories? 
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Recall that anyons are fundamentally defined quantum mechanically and all the energy eigen- 
functions of the system vanish on the set A of coincident points. We modeled the classical system 
by removing A and explored the consequences in the previous sections. In particular we just 
omitted trajectories that could cut through A. 

From a purely classical point of view, though, this is a little unsatisfactory. One could le- 
gitimately ask just what happens to trajectories (e.g. J = for = 2) that attempt reaching 
the disallowed region? To decide this one has to extend the classical modeling by supplementing 
the equations of motions with a specified "boundary" condition. The choice is helped by noting 
that the curl of the vector potential is a 5-function which may be "regulated" to stipulate the 
"boundary condition". This analysis is given in appendix A. 

The result is that all the elliptical trajectories are identical to those of the oscillator and only the 
degenerate trajectories get modified. A degenerate trajectory must reflect back from the coincident 
point. Clearly such a trajectory will have half the period of the generic trajectories and hence we 
do have such "half trajectories". While this can been seen directly in the oscillator confinement 
from the form of the density of states, it may also be seen in the numerical computation of the 
Fourier transform of the density of states of a particle moving on a two-dimensional disc in the 
presence of the flux line|32j. 

While classically there are half-period orbits, the inference that they are manifested in the 
two sub-leading terms of ea.(f75|l. is ambiguous. To see this let us express the propagator in a 
sum-over-orbits form. To do this rewrite the the last two terms to resemble the second and the 
third terms. 

The last two terms also have a (—1)^ in the summation over k and this can be handled in two 
ways. 

(A) Separate these sums into even and odd integer sums. All the sums over k being geometric 
series can be done explicitly {ie needs to be added). Defining 8 = E/{hu)), we get, 



hLog{E) = 8 + 

\{8 + aq-\ )e2^*(^+"^) + ie''^(^+°') 



+ C. C. (78) 



^ 2 y-' 2 1^ if _|_ Q C 

\ — Q2-Kiaq(8-aq) 

Here C. C. means complex conjugation. Now we can "read-ofF' G{E + ie) and get (suppressing 
the Thomas- Fermi term). 



1 1 p2m{€+aq+ie) 

i-^-^)-'^G{E + ^e) = -{8 + aq--{l- e^^^^-^^-)) } {-^^^-^^^} + 

1 1 ■ ^c- , ■ ^ ^2TTi{£-aq+ie) 

2 - + - e-(^--'-^-)) } (79) 
Noting that the second group of braces is sum of a geometric series: 

^2-Ki{£±aq+ie) 



e2.iM£±",+iO, (80) 



1 _ p2-Ki{£±aq+ie) 

^ k>l 

we see that the propagator is expressed in a form suggestive of a sum over periodic orbits. The 
first group within the braces being the "amplitude" while the second group being contributions 
from multiple traversals of a basic periodic orbit. The two terms can be seen to come from families 
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of elliptical orbits going clockwise and anti-clockwise and that there are no term corresponding 
to degenerate ellipses or half orbits. Note that the terms in the "amplitudes" other than f it 
are sub-leading relative to Szizag. This leading part of the "amplitudes" can also be seen to come 
from POT in presence of continuous families of periodic orbits j29| . The "amplitudes" however 
are complex. For = there is cancellation of these sub-leading terms and one recovers the 
exact result for the oscillator. It appears that though the system is integrable, the trace formula 
does not give exact propagator. 

The poles in G{E) come from the two terms and are given by: 

E+{n+) = huj{n+-aq) n+ > 2 (81) 
= huj{n^+ag) n_ > 1. (82) 

The residues at these poles { £±aq = n± ) are n± if n± is even and n± — 1 if n± is odd. These 
residues of course give the degeneracy. Note that the sub-leading terms are important for this. 

It is easy to see that the locations of all the poles can be re-expressed in the form given by 
the exact spectrum. The degeneracies at the poles also match exactly as expected. 

(B) We can also write {—l)^cos{k0) = cos{k{6 ± vr)). Proceeding exactly as in the case (A), 
we can "read-off' G{E + ie) and get, 



\ 27ri ) 



G{E + ie) = + + 

^4 g7ri(£-l+Qq+ie) ^ 4 ^ni(£+l-aq+ie) ^' ^^^^ 



These have poles which are subsumed by the poles from the first two terms. The overall pole 
structure and residues of course are exactly same as before as they should be. 

The exponents in the last two terms however do not look like the action integral. A reflecting 
trajectory will not receive a contribution from the Og dependent total derivative term. The 
reflecting trajectories will also form a single family of trajectories and thus should give only 
one term and not two terms. The last two terms can be combined and the propagator may be 
re-expressed as 



V 27ri J 



G{E + ie) 



E 

k>l 



+ Og ^ ^2mk{£+aq+it) _|_ ^2TTik{£-aq+ie) 



^ ^ -i sin(fc7rag) ^ ^ink{£~i) 



(84) 



The propagator now does have form indicating a contribution from the half period orbits in- 
cluding a contribution from a Maslov index due to reflection. However it has a complex amplitude 
with a dependence on Oc and the multiple traversals index k. 

Thus the exact propagator can be expressed in two different forms mimicking the trace formula 
with and without the half period trajectories. Although the phases are as expected, the amplitudes 
are not. The "amplitudes" are complex in general and also have terms with different orders of h. 

The different orders of h seen can be understood in the context of trace formula as due to 
families of periodic trajectories and not due to any higher order corrections to the SPA. When 
periodic trajectories come in continuous families, the number of Gaussian integrations is reduced 
since the integration over the family must be done exactly l29^ and this increases the powers 
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of h in the denominator. The ehiptic trajectories form a three parameter family while the half 
trajectories form a one parameter family which explains the difference in the powers of h. 

While the powers of h can be understood, the precise matching of the amplitudes from the 
POT does not follow. One may conclude from this that the exact propagator is not of the form 
implied by POT and that additional contributions are necessary. The "half period" peaks in the 
Fourier transform is a property of the exact propagator and need not imply presence of "half 
orbits" in a POT, a "diffractive" contribution may also generate such a peak. 

This analysis of the exact propagator shows that even though the system is integrable, the 
propagator is not expressible as a POT sum, with or without "half orbits" . In the light of the work 
on "diffractive" contributions cited above, the exact propagator seems to need such contributions. 
However we have not done a precise demonstration that the exact propagator can be obtained by 
inclusion of "diffractive" contributions. This is beyond the scope of present review. 

Needless to say, the structure of "half trajectories" will be more complicated for N > 2. One 
needs to develop/adapt the machinery of diffractive contributions to this case to strengthen a 
semiclassical demonstration of the non-linearly interpolating eigenvalues. 

Observe that even without the knowledge of precise number of families, we have precisely 
2N{N-i)/2 g^j^g ^}^g sum over periodic trajectories. This is because, the contribution of the 
action integral depends only on the orientation class and not on further possible subclasses of 
trajectories. The measure factors therefore add up within each orientation class and this is 
sufficient to get the semiclassical eigenvalues. 

The sum over families in the trace formula contains as many terms. To get the eigenvalues 
we need to look at only the action integral which differs from that of the oscillator only by the 
additional <f ^ Yli^j ^ij contribution. The oscillator's contribution to the action integral is of 
course just "^E, where E is the classical energy and uj is the oscillator frequency. 

The net Oc dependent contribution to the action integral is given by. 



Here eij is the sign of traversal of fij and is ±1. The is trivial to evaluate. 

Since the classical trajectories are identical to those of the oscillator, we will have the same 
Maslov indices which incorporate the usual zero point energy namely, A^(A^— 1) for 2N dimensional 
oscillator. Including this contribution, the semiclassical energies are given by : 



The coefficient of takes all possible integer values from —A^(A^—l)/2 to N{N—l)/2 insteps 
of 2. The extreme values correspond to all the particles traversing the same way and represent the 
contribution of "collective" motion shown in The remaining values represent contributions 
of "relative" motion. These are new eigenvalues albeit at the leading h level. Numerical spectra 
available for = 3, 4 show eigenvalues matching with the above at the = and Uq = \ 
|22j . These numerical eigenvalues however are non-linearly interpolating while our semiclassical 
eigenvalues are linearly interpolating. 

In reference |22j also it was conjectured that there will be linearly interpolating eigenvalues 
with various slopes based on a different semiclassical argument. We have obtained these eigen- 
values by direct application of the periodic orbit theory and proved their existence for all A^. 




(85) 





n > 0. 



(86) 
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Observe that the Og dependence in the eigenvalues ( locations of poles of G{E) ) arises only 
from the explicit ac ^ cxqh in the action. This therefore must be at the most linear in ag. Since Og 
is fixed, the semiclassical limit of ^ — implies ac ^ 0. Thus one may at the most see linear Oc 
dependence in the leading approximation which could however be indicative of exact non-linearly 
interpolating eigenvalues. 

Thus the leading level semiclassical propagator (with POT only) already indicates the presence 
of eigenvalues such that E{ag = 1) — E{ag = 0) takes all possible integer values between ±N{N — 
l)/2 in steps of 2. 

As there does not appear to be any scope for going beyond the leading SPA (all intermediate 
exponents of phases are quadratic), it may be conjectured that diffractive contributions are needed 
to improve these semiclassical eigenvalues. 

5.3 Pseudo-integrability 

In the previous sections we encountered two important features, one related to the fundamental 
group and one related to symmetry reduction due to incomplete vector field. Both were caused 
by the same source, namely removal of coincident points implying topologically non-trivial phase 
space. Both are relevant for periodic orbit theory in a general way and a few general remarks are 
in order. 

The stationary phase approximation to the propagator naturally leads to periodic orbits in 
the phase space. These orbits may be isolated or come in continuous families or both. The con- 
tinuous families may be generated by a full group of symmetries or by only a subset of symmetry 
transformations. Since each periodic orbit is also a map of to the phase space T, clearly every 
orbit must belong to one and only one homotopy class of the fundamental group, vri(r), of the 
phase space. If an orbit is a member of a continuous family then the entire family must belong 
to a single homotopy class. Note that a given homotopy class may contain no orbit (solution 
of equation of motion) or several isolated orbits and/or several families of orbits. But a family 
can not spill over two distinct homotopy classes. In sub-section 5.1 above, we saw precisely the 
splitting of a single basic family for oscillator into two basic families for anyons because of the 
non-trivial fundamental group. Multiple traversals of course belong to different homotopy classes 
and are explicitly summed over. Thus a non-trivial vri may (but not necessarily) provide an ob- 
struction to a symmetry. How exactly may such an obstruction manifest itself? For this we have 
to consider vector fields generating symmetries. 

Recall ^33' that every function on the phase space generates infinitesimal symplectic diffeomor- 
phisms (canonical transformations) via its corresponding (globally) Hamiltonian vector field. The 
Lie algebra of such vector fields is isomorphic to the Poisson Bracket (PB) algebra of functions 
on r. However such infinitesimal transformations exponentiate to give a one parameter group of 
transformations only if the vector field is complete i.e. the integral curves of the vector field can 
be extended so as to have these as a map from the full R. Only complete vector fields — and this 
is a global statement — give rise to groups of symmetries. (A corresponding quantum mechanical 
statement for continuous symmetries is contained in the Stone's theorem: every one parameter 
group of unitary transformations is generated by a self-adjoint operator and conversely.) Since we 
usually want classical symmetries to be reflected at the quantum level with observable generators 
we have to have groups of symmetries and hence complete vector fields. 

The criterion of integrability in terms of vanishing PB's guarantees only the existence of 
infinitesimal symmetries which is of course a prerequisite. An "infinitesimally integrable" system 
may thus be: (a) integrable via action-angle coordinates if the vector fields are complete and the 
integral curves are closed; (b) integrable, but not by action- angle variables if the vector fields are 



27 



complete but only a subset of these have closed integral curves and (c) partially integrable if only 
a subset of vector fields are complete. The "pseudo integrability" property of anyons pointed out 
in ref. ^H] falls in the category (c). As explained in the appendix B, apart from the Hamiltonian 
only the total angular momentum has a complete vector field and hence is the only quantum 
number that survives. 

A non-trivial fundamental group by itself however does not imply possibility of incomplete 
vector fields. For, on a compact manifold all non-singular vector fields are complete j33j and it 
can of course have a non-trivial tti. In many cases constant energy surfaces are compact and one 
does not have to worry about incomplete vector fields. For the anyons however removal of the set 
A makes the constant energy surface non-compact which admits possibility of incomplete vector 
fields and corresponding loss of symmetry. In the above, we have seen manifestation of all these 
features. 

6 Semi-classical analysis of the ground state 

The semiclassical analysis of the last section is a general one for all the states of the N-anyon 
system. However, the ground state of this system deserves special attention, being somewhat 
unusual. As pointed out by Chitra and Sen [2^, due to the level-crossing between pairs of energy 
levels, the ground state consists of many pieces at different a values. This is because states of 
different J have lowest energy at different values of alpha. In this section we present briefly a 
semiclassical analysis of the many-anyon ground state via another approach, which is motivated 
by the above observation. 

We start from the classical solutions of lowest energy, which are the stable configurations 
("classical ground states") of the system. These states are determined by minimizing the total 
energy of the system. The assumption is then made that in the semiclassical regime, the quantum 
ground state energy is given by the classical energy to the zeroth approximation with corrections 
to first order given by quadratic, quantum fluctuations about this value. 

Now the quantum ground state doesn't interpolate smoothly between the extreme values of 
the statistical parameter. This behaviour can never be captured if we try to find the ground state 
by simply minimizing the classical Hamiltonian. In fact, such a step gives a "ground state" which 
is completely independent of the angular momentum, and therefore of a. 

The rotational symmetry of the Hamiltonian also means that the angular momentum of the 
system is conserved. To find the minima of the energy, one can proceed in two steps: first 
minimize the Hamiltonian for a fixed angular momentum and then minimize the resulting energy 
with respect to the angular momentum. In the second step, one could carry out the minimization 
over quantized values of the angular momentum. This method, as we shall shortly see, retains the 
dependence on angular momentum of the lowest energy, even at the classical level and also has 
some features observed in the Thomas- Fermi approach to the ground state |21j . 

6.1 Classical ground states 

The energy of N free anyons in a harmonic confinement is represented by the classical hamiltonian 




where 




(87) 
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Constancy of angular momentum is applied as a constraint, 



N 



C = ^ji - J where ji = XiPy^ - yip^ 



Introducing the Lagrange multiplier A, the constrained function 

F = H -\C, 

must be minimised. That is, 5F = gives the extrema which will be (local) minima if > 0. 
The equations for equilibrium are thence derived to be the following: 

Px, = [ax, - A)2/i, 

Py^ = (ay^ + A)xi, (90) 

are conditions arising from the variation of F with respect to the momenta, while variation with 
respect to the positions gives the equation 

(1 - X^)r, = 0. (91) 

Since this is true for all i, this fixes the Lagrange multiplier: A = ±1. The constant angular 
momentum becomes 

N 

J = -X^rf-^N{N-1). (92) 

i 

The equilibrium configurations are thus dynamic paths with non-zero momentum and angular 
momentum, the positions being given by those satisfying the above relation for a given angular 
momentum sector. The energy in such a sector is given by 



N 



= ±{J + ^N{N-1)). (93) 

To find the ground state energy £"0, one must now minimize this energy with respect to J. Since 
the main aim of the present exercise is to find a semiclassical route to anyon ground state energy, 
the above expression can be minimised with respect to the quantized, integer values of J. Now 
the term ^N{N — 1) can be expressed as a fraction plus an integer: fj, + n where n G Z. So the 
energy is quantized as Eq = 171 + fj, where m is an integer and /i is a fraction. Eq is evaluated and 
plotted as a function of a for different N in figure ^ 

A remarkable feature of this picture is that the the ground state for different values of a has a 
different J. Note that the ground state energy for = 2 is exact, even quantum-mechanically. For 
higher A^, the "sawtooth" envelope of the minimum energy over the whole range of a reflects the 
behaviour observed in the earlier numerical studies on the quantum ground state energy ^1 . 
Level crossing effects are seen at the classical level itself, as different J values contribute to the 
ground state at different a. Quantum corrections could be expected to smooth out this graph. 
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Figure 1: Classical ground state energy of two and three anyons as a function of a 
6.2 Towards Semiclassical corrections 

Quantum mechanically, the ground state energy will be corrected by the zero-point energy due to 
fluctuations about the classical minimum energy configuration. Now we can make the assumption 
that the potential is approximately harmonic near the minimum positions. The quantum wave 
function is localized around these points. To calculate the corrections to the energy, one can 
expand the potential as a Taylor series about the minimum. The zeroth order term is of course 
the classical ground state energy. The first order term is zero and the second order term is 
quadratic in the fluctuations. The matrix of second derivatives of the potential, or the Hessian, 
can be diagonalised in the basis of the normal modes of the system. Then the second order 
term for the expanded potential is just a set of oscillators in the normal modes, whose squared 
frequencies are the eigenvalues u>f of the Hessian. These oscillators can be quantized to give 
corrections AE = 1/2 '}2iLi to first order in h. This of course is valid only if one can neglect 
the higher order terms in the expansion of V (see, for example, Raiaraman|34j. chapter 5). One 
hitch in this program is the occurrence of zero modes: when any one of the Wj's is zero, then this 
approximation is no longer valid. This happens when the configuration is one of neutral stability, 
i.e. there exist directions of symmetry in the model. Perturbations along the symmetry direction 
do not change the energy of the system. The quantum energy eigenfunctions cannot be said to be 
localized near the minimum: they will spread along the flat "valley" . Zero modes can be dealt with 
in various ways. One of the most straightforward methods is discussed in chapter 8 of |31|. The 
basic idea is to work in a system of coordinates where the symmetry directions get separated when 
calculating the quantum corrections. The method we will adopt is to perform a transformation 
in the coordinates such that the direction of symmetry as the collective coordinate becomes 
apparent and can be isolated. The Hamiltonian is manifestly independent of this coordinate 
and the corresponding conjugate momentum is conserved. Hence the fluctuations along this 
direction can be ignored and the Hessian, being only in terms of the remaining coordinates, picks 
up the contribution from fluctuations along directions that are not unchanged by the symmetry 
transformation. This is equivalent to transforming the Hessian to a basis in which it looks block- 
diagonal, such that one block contributes the zero eigenvalues alone and the others do not. This 
block can be separated out. One then diagonalises the remaining part of the matrix. 
In the N-anyon case, we work with the Lagrangian 




(94) 
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where is the usual statistical gauge potential and V(ri) is some confining potential. 

Overall rotations is a symmetry and motion along the rotation coordinate can be treated as 
a collective coordinate and separated out. This was already seen in section 4: we perform the 
time-dependent transformation 

= R{9)rji with 7]y^ = 0, (95) 

The Lagrangian then becomes 

N N N 

1 11 
The conjugate momenta are defined by the equations 

Px, = r]x, - 0r]y^ + (98) 
Pyi = rivi + Or],^ + Ay„ (99) 
PxN=rixN+AxN^ (100) 

N N-1 

i=i 1 ^ 

It is convenient to perform a canonical transformation to new canonical momenta 

ni = Pi- Ai. (102) 
Taking = X, the Hamiltonian is then 

N / N-1 \2 

^ = 2 + + f ^(^ - 1) - 77, X ttJ + Vin). (103) 



and constrained minimization gives 



TTx = 0, 



(104) 



For harmonic confinement, V = ^J2Vi ^^id the Lagrange multiplier gets fixed to ±1, the choice 
depending on which minimizes the energy. The equilibrium conditions can then be written down 
as 

TTi =Ak X rji, 

i 

This means that the classical extrema lie on a (2A^ — l)-dimensional hyper-sphere in coordinate 
space whose radius is the square root of the magnitude of the total angular momentum. The 
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energy is 



N 

2 



=\J\ ^^^^^ 
=\P, + ^N{N-1)\ 

which is the same as obtained before. The advantage is that in these coordinates the Hessian 
will not contain the zero mode corresponding to rotations. The Hessian matrix, evaluated at the 
equilibrium positions, in terms of the coordinates = 1 . . . A'' — 1, momenta TTxiiT^yi^i = 

1 . . . N — 1 and finally the coordinate X and momentum ttx of the Nth. particle, looks like 



M = 2C^ 4 0, (107) 




where 



A 



X 

I|2iV-2 ^\2N-2 
-^\2N~2 I|2Af-2 



Q = ( 1^ the symplectic form in 2N — 2 dimensions. 

Y-i|jV-i U / 

This matrix must now be diagonalised to get a quadratic form in the normal modes. In this 
case, since the ground states arc not static this matrix is in terms of all the phase-space variables 
and not merely the coordinates. The diagonalization must therefore preserve the symplectic 
structure of the phase space. The details of this symplectic diagonalization are given in appendix 
C. 

We find that the matrix has one eigenvalue equal to 1, 2N — 2 eigenvalues equal to 2, one eigen- 
value equal to 2 ( 1 + ) and finally, extra 2N — 2 zero modes. The eigenvectors corresponding to 



these extra zero modes are tangential to the equilibrium phase space coordinates. Displacement 
along these directions preserves the 2N — 1 dimensional sphere of equilibrium, so they correspond 
to displacements along this sphere. These zero modes cannot however, be removed using similar 
techniques as used above. The main reason for this is that they correspond to infinitesimal dis- 
placements along the 2N — 1 dimensional sphere mentioned above, which is however punctured at 
the positions of coincidence of particles. The displacements cannot be made large, and therefore 
these zero modes do not actually represent a global symmetry direction of the system. 

This path to finding the exact ground states semiclassically seems to face difficulties here, and 
the origin seems to be precisely the analytic intractability of the non-trivial energy levels that will 
descend through level crossings. 



7 Summary and Discussion 

We have presented in this review several issues related to the classical and quantum mechanics 
of anyons. Since several books and reviews dealing with various aspects of anyon physics already 
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exist, we have focussed on aspects that have received httle or no attention until now. With 
this in view, we have specifically focussed on the integrability properties of the anyon system 
and semiclassical analysis of the spectrum. We have therefore discussed two properties of the 
many anyon system: (1) The partial separability of collective and internal degrees of freedom 
and (2) its identification as a pseudo-integrable system. We have shown that the first property 
explains the existence of the exactly known eigenvalues which is some what uncommon for a 
generic many body system. This also shows that the exactly known spectrum incorporates only a 
some what trivial aspect of anyon dynamics. The nontrivial aspects although partially uncovered 
by numerical results are still remain elusive. 

The second property enables one to understand qualitatively the origin of regular and irregular 
features in the spectrum of many anyons. We have clarified and sharpened the meaning of pseudo- 
integrability. The system is locally identical to the isotropic oscillator in two dimensions but not 
at the global topological level. It is the non-trivial global topology that reduces the symmetry 
group from that of the oscillator. 

We also considered application of the periodic orbit theory to this locally trivial (integrable) 
but globally non-trivial system. As a by-product, we saw how the exact spectrum for two anyons 
is reproduced by semiclassical methods. We saw that while the dynamical symmetry is reduced 
from SU{2) to U{1), the rank remained the same and hence integrability property is preserved. 

For > 3, we reproduced the previously known exact eigenvalues JHI- In addition, we 
obtained further new linearly interpolating eigenvalues. In the language of reference these 
correspond to the signature of 'relative' dynamics. In this case the symmetry reduction was 
drastic, from 0Sp{4N,R) to OSp{i,R) x 0{2,R) x 0(2, R). The rank was reduced from 2N to 
6, destroying the integrability property. 

We have also discussed in detail the issue of "half trajectories and pointed out their ambiguous 
role. Taking the view-point that "half trajectories be excluded, and noting that there does not 
seem to be any scope for computing higher order corrections to the semiclassical spectrum with 
a possible non- linear dependence on aq, it seems that the non- linearly interpolating eigenvalues 
are genuine quantum consequences beyond what semiclassical analysis could give. Semi-classical 
analysis is nevertheless sufficient to indicate the presence of these eigenvalues. Thus, the impact of 
pseudo-integrability of the anyon system on the classical periodic orbit theory is quite non-trivial. 
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Appendix A : Regularised classical dynamics and reflecting orbits 



Consider without loss of generality the case of two anyon in the relative coordinates. Generic 
orbits in the configure space, B? — {0}, are of course an ellipse and only a degenerate elliptical 
orbit attempts to pass through the origin. Noting that this system can be thought of as a charged 
particle in presence of a singular ("statistical") magnetic field at the origin, one may regularise 
the magnetic field or the flux to study the orbits and obtain the limiting behaviour to deduce 
"boundary condition" at the origin. 

To do this we imagine the relative Hamiltonian to be that of a particle in a magnetic field 
along the z axis. The field is of course to be effectively confined to a small disc around the origin. 
Now observe that for an axially symmetric magnetic field along the z-axis, -B(r, 9) = B{r)k where 
r, 9 are the usual spherical coordinates in two dimensions, we may write the vector potential as, 

A = Arr + Ae9 (108) 

This implies, 

B{r) = drAe + — - -deAr, (109) 
r r 

In a symmetric gauge, the vector potential is independent of the polar angle and by choosing 
it to be divergence free one can set the radial component to zero. 
The fiux ^{R), through a disc of radius R is given by, 

$(i?) = 2vr^ clrrBir) = R A0d9 = 2TTRAg{R). (110) 

This implies, 

Ae{r) = ^ = -l dr'r'B{r') (111) 

The choice, <I>(r) = 2iTa ^ r > 0, gives the vector potential used in the quantum mechanical 
calculation. It also implies the rB{r) = a5{r) and this of course is the singular nature of the 
magnetic field. Notice that smearing the 5{r) will not make the magnetic field non-singular at 
the origin because of the explicit 1/r. 

For a regulated system one wants the magnetic field to be non-singular every where and 
effectively confined to a disk of radius e around the origin. This requires the vector potential Aq 
also to be non-singular and therefore vanishing at the origin. If Ag{r) cr^ as r approaches 
zero then, B{r) c{/3 + l)r^^^. For a finite, nonzero B{0) we must have /3 = 1 and therefore 
$(r) 27r cr^ near the origin. For r > e we still retain the flux to be 27ra. Continuity at r = e 
then gives c = a/e^. Notice that this limiting behaviour is fixed by the demand of non-singularity 
of the fields and as such must be reflected in any explicit choice for the magnetic fleld. 

Thus our regulation involves choosing, 

$(r) = 27ra V r > e 

^ 27ra^ as r ^0 (112) 

One could choose a uniform nonzero magnetic fleld inside the disk as an explicit choice but 
it will not be necessary. Since we are interested in the limiting behaviour of trajectories as e is 
taken to zero, the limiting behaviour of the flux is all that we need. 
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Consider now the orbit equation for e nonzero. Integrating the equations of motion once using 
the two constants, energy E and angular momentum £, the orbit equations in r, 6 coordinates 
become : 



r 



r^e = £-rAg{r) (113) 




r < e 
r > € 



There are three types of orbits possible: those which are fully inside the disc, those which are 
fully outside the disc and those which go both inside and outside the disc. In the limit of e going 
to zero, the first type of orbits are clearly irrelevant. It is easy to see that for the second type of 
orbits one must have i a. These are insensitive to the flux in the limit and are thus identical 
to the orbits of the oscillator. 

For the last type, an interior turning point is possible only for < i < a. The condition 
that such an orbit must also have an exterior turning point limits i to a. We are interested in 
computing the change in the angular coordinate from the entry into the disk till exit from it. 
Explicit computation shows that the change in the angular coordinate goes to zero as e goes to 
zero. Thus such an orbit reflects at the origin. Note that these are precisely the radial orbits. 

To summarize, a regulated classical modeling for anyons is generically stipulated by giving the 
behaviour of the flux near the origin. It amounts to cutting out a disk of radius e and filling it 
up with non-singular fields. The classical non-radial orbits then are exactly same as those of the 
oscillator except for the replacement I I — a. The radial orbits reflect at the the origin and are 
termed half orbits since their period is half of that for the other orbits. 
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Appendix B : OSp(4N, R) dynamical symmetry and classification 
of trajectories 

In this appendix we collect together some of the well known facts about the dynamical symmetry 
of the n dimensional isotropic oscillator. 

The isotropic oscillator in n dimensions has OSp{2n, R) as the group of dynamical symmetries. 
This is a group of 2n x 2n order matrices, g, which are both orthogonal and symplectic. Denoting 
by g the transpose of g, we have the defining equations: 

g g = hn ■ (orthogonal) 

g g = ri : (symplectic) 
and with g + eT generators T satisfy 

T = -T 

T n = -nT 

Here is the identity matrix of order 2n while f2 is a suitable matrix defining the symplectic 
condition. This will be chosen below. 

The symplectic condition ensures that wc have a (linear) canonical transformation while or- 
thogonality ensures that the Hamiltonian, J2iiPi + , is invariant. It is easy to see that 
dimension of this group is while its rank is n. In fact one can show that the group OSp{2n, R) 
is isomorphic to the group U{n, C). Making an explicit choice of in a block form, one can obtain 
the block form for g using the symplectic condition. The real block matrices can be combined into 
complex block matrices. The orthogonality condition in terms of real matrices then translates 
into the unitarity condition for the complex matrices thereby proving the group isomorphism. 
This isomorphism immediately implies that the ortho-symplectic group acts transitively on the 
constant energy sphere. This is used in the section 5.1 . 

We are interested in finding the action of this group on the phase space via canonical trans- 
formations (symplectic diffeomorphisms) . For notational convenience let us group the standard 
canonical variables together and denote them by cj^ , /i = 1, 2, 2n, the first n being coordinates 
and the last n being the momenta. The dual of the symplectic form has components, O'^'^ given 
by the iiu th element of the matrix fl. With this notation, the PB of functions on the phase space 
and the infinitesimal canonical transformations are given by, 

{Fiio) , G{io)} = n>"' d^F 

S.uj^' = e n^" d^F{uj) (114) 

Functions purely quadratic in to generate linear canonical transformations and are also closed 
under the PB's and those which leave the Hamiltonian invariant give the symmetry transforma- 
tions. Explicitly, 

F = ^Cd A u , G = ^io B Lo ^ 

{F, G} = \u: (AfJB - BfiA) 

The Hamiltonian corresponds to A = l2n- The Poisson bracket of any G with B. vanishes 
provided the matrix B commutes with 17. Thus generic matrices defining quadratic functions 
are real, symmetric matrices commuting with f2. The matrices f2B are then antisymmetric and 
provide an isomorphism of quadratic functions to the generators T of the group 0Sp{2n, R). 

Action of the one parameter group generated by a function G is found from the integral curves 
of the corresponding Hamiltonian vector field. These curves are defined by the matrix equations, 
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(T = QB) 



^ = QBu (115) 
da 

u{a) = (e'^"^)a;(0) (116) 

It is convenient to choose a particular grouping of the phase space coordinates and correspond- 
ing choice of ft. Firstly let us put n = 2N, relevant for the present context, so that the phase space 
is 4:N dimensional. Arrange the coordinates and momenta as xi,yi,pix,Piy, ■■■■,3Sn ,yN,PNx,PNy 
and denote by coi the coordinates and momenta of the i^^ particle. The index i now runs from 
1, N and each cjj is a 4 x 1 matrix. Correspondingly we choose f2 as a block diagonal matrix 
with N blocks and each block being 4x4 matrix A. We choose, 

X ll ) 

One can choose a basis for OSp{AN, R) as follows. Let us denote the generators as Tj, 
i = 1,2, ...,N and Ty with i < j. Each of these are expressed in the block form. The Tj are 
block diagonal with a nonzero 4x4 matrix Uj as the i*^ block element. The Tij have a matrix 
Vij at the z*'* row and j*'* column {{ijY^ block) and — Vy at the (ji)*'^ block. In equations. 



With these definitions it is easy to translate the conditions on the generators in terms of the 
4 matrices u , v as: 

Uj = - Uj UijA = Auij ; A = Avij (119) 

Thus the u's generate OSp{A, R) while the v's are required to commute with A. The number 
of independent u's is 4 while number of independent v's is 8 of which 4 are 'diagonal' and 4 
are 'off-diagonal'. The dimension of OSp{AN,R) is thus AN + 8N{N - l)/2 = 4N^. These 
independent matrices can be explicitly chosen in terms of 2 x 2 Pauli matrices and the identity 
matrix as: 

/ ia2 O2 \ f O2 h\ f O2 ai\ f O2 as 
~ ^ O2 .a J ' V -I2 J ' V -^1 J ' V -^3 O2 
f h 02\ f ai 02\ f ia2 O2 \ ( CJ3 O2 \ 
""^''"''''^ ~ [02 I2) ^ [02 aj ^ [02 ^a2) ' [02 as) ' ^'''^ 

O2 I2 \ ( O2 ai\ [02 ia2 \ ( O2 a3 
V(5,6,7,8) ~ \-l2 O2 J ' \ -ai O2 J ' \ -ia2 O2 J ' \ -as O2 

We are interested in one parameter groups generated by some element T of the Lie algebra. 
If T satisfies: = -P with P satisfying = P , [T, P] = 0, then it follows that, 

e""^ = I4JV - P + P [cos (a) sin (o-)T] P (121) 

For the basis generators Tj we have P to be block diagonal with I4 as the i*'' block while for 
Tij we have P to be the block diagonal matrix with I4 as the i*'* and the j*'* blocks. Using these. 
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the exponentials of basis generators can be evaluated to get the integral curves as, (with obvious 
notation) 



uj{a) = [ I ■ + Ij { cos{a) + sin{a) Tj } Ij ] uj{0) and, 

uia) = [ I'ij + lij { cos{a) + sin{a) } lij ] lo{0) (122) 

Therefore integral curves of the basis generators are periodic curves. Further, the Tj's affect 
only the z*'* particle position and momenta while the Tjj mix the z*'* and j*'* particles only. 

We need to study how the angular momenta Ji and Jij vary under the action of one parameter 
subgroups. Define the (block) matrices, 

(Lj)mn = ^SimSmn L = V7 

In terms of these matrices the angular momenta are given by: 

Ji = ^LoLi LO = ^UiL Ui 

Under the group generated by the basis generator, Tj, the J^ for instance varies as: 

2Jm(o') = 2J„ + Sim [ -sin^^a) {2Jj - 0)1 (u L u) ui} 

+ sin{a)cos{a) { u)\ (uL + Lu) oj^ } ] (125) 

while under the group generated by the basis generator, T.^, the Jm varies as: 

2Jm(o") = 2J^ + 5im [ -sin^ia) {2Ji - uij (v L v) wj} 

+ 2sin{a)cos{a) { lD; (Lv) } ] 
^jm [ —sin^{a) {2Jj — (v L v) uji} 

— 2sin{a)cos{a) { cDj (vL) } ] (126) 

Similar but more complicated expressions follow for the Jmnio') also. 

Remark: For the case of = 1 we have only one block. This could be either a single two 
dimensional oscillator or the relative coordinate dynamics of two anyons. The generators are 
of course only u's. Of these U2 corresponds to the Hamiltonian itself while Ui generates same 
rotations of both r, p. These two matrices commute with L while the remaining two anti-commute 
with L. This leads to the result quoted in the section 5.1 . 

To deduce the surviving symmetry group for many anyons we need infinitesimal variations of 
Jmn for arbitrary generator T. This is easily derived and is given by, 

SJmn = ^tl^mn L (T^^- - T^j )u;j (127) 

j 

In section 5.1 we needed to determine T such that SJmn is zero whenever Jmn is zero. That is, 
at all points in the phase space where any single Jmn is zero, we want its infinitesimal variation 
induced by T to be zero. 

Fix a particular m and n. Fix tOmn- We can consider points with Uj, j ^ m,n such that 
no other Jij is zero. SJmn = then implies Tmj = T^nj V j 7^ m, n. We can also consider 
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points with different uJm,uJn keeping LOmn fixed and maintaining all other conditions. This implies 
Tmm — T^nm + T„j„ — T„„ = 0. Repeating this for all m ^ ra fixes the form of T in terms of 
arbitrary generators of 05p(4, i?), u , v as, 

(T)y = uSij + v(l - (5y ) , [ L , u - V ] = (128) 

The choice u = v corresponds to Tjj = u V i,j. It follows that all Jmn 's are invariant 
independent of their values. It is easy to see that these four T's effect transformations of the 
center-of-mass variables which are insensitive to the anyonic features. This OSp{A, R) symmetry 
is thus always present for all A*" > 2. 

The choice v = implies T is block diagonal with the same u on all the diagonal blocks. 
Further, [L, u] = implies that u must be a linear combination of ui and U2 . These two 
T 's can be seen to correspond to the total Hamiltonian and the total angular momentum. 
This result is used in the section 5.1 to deduce that the surviving symmetry for N anyons is 
OSp{A,R) X 0(2, i?) X 0{2,R). 
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Appendix C : Symplectic Diagonalization 



In our method of calculating semiclassical corrections to a classical ground state, we need to 
diagonalise the Hessian matrix, the matrix of second derivatives of the potential at the ground 
state.. In the case of anyons, since the ground state is not static this matrix is in terms of all the 
phase-space variables and not merely the coordinates. The diagonalization must hence preserve 
the symplectic structure of the phase space. The Hamiltonian expanded about the ground state 
is 

H = Hq + (^Mi + higher order terms, (129) 

where ^ is the column vector made up of the phase space variables. 

The phase space structure is preserved if the matrix M is diagonalised by a symplectic matrix 
S such that S'^IS = J, J being the 2N — 2-dimensional symplectic form. Now there exists a 
theorem^ known as Williamson's theorem 36 (a compact summary is presented in Appendix 6 of 
Arnold |221) part of which implies that if M is a real symmetric positive definite matrix then it can 
be symplectically diagonalised by such an S. Moreover, if S'^MS = where D = diag(Di, D2), 
Di,D2 being — 1 dimensional diagonal matrices, and M = — tM~^/^JM~^/^ is a Hermitian 
matrix with real eigenvalues ±uji then it can be shown that D2 = diag(wj). When this 
happens, then if the normal modes are ^' = ( p' ) = S~^^, the quadratic form looks like 

i (130) 

This is a harmonic oscillator in the space of normal modes, on quantizing which we get the 
correction to the ground state energy as 

l^E.-^r'- (131) 

Thus one needs just the positive eigenvalues of the matrix M, which should further be non-zero. 
This will require M itself to have non-zero eigenvalues. 

Now to find the eigenvalues of the Hessian (|107|) : The last row and column, coming from 
the derivatives with respect to Px, contribute an eigenvalue 1. Consider the reduced 4A^ — 3 
dimensional matrix M' formed by removing the last row and the last column. It turns out that 
this matrix does have zero eigenvalues, inherited from the eigenvalues of A with eigenvectors of 
the form k{c, -^C^c) , k being an arbitrary proportionality factor. There are exactly 2A^ — 2 
of these. The remaining eigenvalues consist of an eigenvalue equal to 2 which is 2A^ — 2 fold 
degenerate and a non-degenerate eigenvalue equal to 2(1-1--^). 



^We are indebted to Prof. R. Simon for bringing this theorem to our notice, as well as for the simple proof in 
the paper 
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